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Abstract. We prove that the Hartree-Fock orbitals of pseudorel- 
ativistic atoms, that is, atoms where the kinetic energy of the elec- 
trons is given by the pseudorelativistic operator \/ ~ A + I — 1 , are 
real analytic away from the origin. As a consequence, the quan- 
tum mechanical ground state of such atoms is never a Hartree-Fock 
state. 

Our proof is inspired by the classical proof of analyticity by 
nested balls of Morrey and Nirenberg 27 . However, the technique 
has to be adapted to take care of the non-local pseudodifferential 
operator, the singularity of the potential at the origin, and the 
non-linear terms in the equation. 



1. Introduction and results 

In a recent paper [5] , three of the present authors studied the Hartree- 
Fock model for pseudorelativistic atoms, and proved the existence of 
Hartree-Fock minimizers. Furthermore, they proved that the corre- 
sponding Hartree-Fock orbitals (solutions to the associated Euler-La- 
grange equation) are smooth away from the nucleus, and that they 
decay exponentially. In this paper we prove that all of these orbitals 
are, in fact, real analytic away from the origin. Apart from intrinsic 
mathematical interest, analyticity of solutions has important conse- 
quences. For example, in the non-relativistic case, the analyticity of 
the orbitals was used in [131 l2l] to prove that the quantum mechan- 
ical ground state is never a Hartree-Fock state (or, more generally, 
is never a finite linear combination of Slater determinants). A direct 
consequence of our main regularity result is that this also holds in the 
pseudorelativistic case. Our proof also shows that any /7^/^-solution 
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<y9 : — )■ C to the non-linear equation 

(V - A + - ± * I • |-i)<^ = Av? (1) 

which is smooth away from x = 0, is in fact real analytic there. As will 
be clear from the proof, our method yields the same result for solutions 
to equations of the form 

(-A + m)V + \/(^+|(^|V = A(^, (2) 

where V has a finite number of point singularities (but is analytic 
elsewhere), under certain conditions on m, s, V, and k (see Remark 11.21 
below). We believe this result is of independent interest, but stick 
concretely to the case of pseudorelativistic Hartree-Fock orbitals, since 
this was the original motivation for the present work. 

We consider a model for an atom with N electrons and nuclear charge 
Z (fixed at the origin), where the kinetic energy of the electrons is 
described by the expression y^(jp|c)M-~(wM5^— mc^. This model takes 
into account some (kinematic) relativistic effects; in units where h = 
e = m = 1, the Hamiltonian becomes 

N 

j=l l<i<j<N ' 

with T(p) = E{p) - a-^ = v^p^To^ - a'^ and 1/(x) = Za/|x|. 
Here, a is Sommerfeld's fine structure constant; physically, a ~ 1/137. 

The operator H acts on a dense subspace of the iV-particle Hilbert 
space Tip = A^]^L^(]R^) of antisymmetric functions. (We will not con- 
sider spin since it is irrelevant for our discussion.) It is bounded from 
below on this subspace if and only if Za < 2/7r (see ; for a number 
of other works on this operator, seePElElIISlISSlESlEIlES). 

The (quantum) ground state energy is the infimum of the quadratic 
form q defined by H, over the subset of elements of norm 1 of the cor- 
responding form domain. Hence, it coincides with the infimum of the 
spectrum of H considered as an operator acting in Tip. A correspond- 
ing minimizer is called a (quantum) ground state of H. 

In the Hartree-Fock approximation, instead of minimizing the qua- 
dratic form q in the entire A^-particle space Hp, one restricts to wave- 
functions "if which are pure wedge products, also called Slater determi- 
nants: 

^(xi, . . . , xjv) = -L det(M,(x,-))5=i , (4) 
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with {ui}fL^ orthonormal in L^(R'^) (called orbitals). Notice that this 
way, e Tip and ||\E'||2,2(i83jv) = 1. 

The Hartree-Fock ground state energy is the infimum of the quadratic 
form q defined by H over such Slater determinants: 

E^^{N, Z, a) := inf{ q(*, | * Slater determinant } . (5) 

Inserting \E' of the form in (j4]) into q formally yields 



In fact, Ui G H^^'^(R^), 1 < i < N, is needed for this to be well-defined 
(see Section [3] for a detailed discussion), and so (I5])-(I6]) can be written 



E"^(iV, Z, a) = inf { £«^(mi, ...,un)\ (mi, . . . , mjv) G Mn} , (7) 
Mn = { {uu . . . , Uiv) G [H'/\M.^)f I {u,, u,) = 6,, } . (8) 



Here, ( , ) denotes the scalar product in L^(M^). The existence of 
minimizers for the problem ([7])-([8]) was proved in ^ when Z > N ~ 1 
and Za < 2/7r. (Note that such minimizers are generally not unique 
since S^^ is not convex; see [10] )• The existence of infinitely many 
distinct critical points of the functional S^^ onAiN was proved recently 
(under the same conditions) in [7]. 

The Euler-Lagrange equations of the problem ([7])-([HD are the Har- 
tree-Fock equations, 



Here, the £j's are the Lagrange multipliers of the orthonormality con- 
straints in dH]). (Note that the naive Euler-Lagrange equations are 
more complicated than ([9]), but can be transformed to ([9]); see [TO].) 






(9) 




N 
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Note that ([9]) can be re-formulated as 

h^pifi = Eiifi , l<i<N, (10) 

with hip the Hartree-Fock operator associated to (p = {(fi,...,ipN}, 
formally given by 

hipU = [T(— iV) — V]u + aR^pU — aK^pU , (11) 

where R^pU is the direct interaction^ given by the multiplication opera- 
tor defined by 

^ Jr3 |x - y| 

and KipU is the exchange term, given by the integral operator 

(/f,„,)(x) = f:(/ «f)dy)^,(x). (13) 

j^-^ ^ jR^ i-^ ~ y\ ' 

The equations ^ (or equivalently f lTU]) ) are called the self- consistent 
Hartree-Fock equations. One has that (Jess{hip) = [0, oo) and that, when 
in addition N < Z, the operator h^p has infinitely many eigenvalues 
in [— a~^,0) (see ^ Lemma 2]; the argument given there holds for 
any = {ipi, . . . ,ipn}, (pi e if^/^(M^), as long as Za < 2/7r). If 
{ipi, ... ,ipiy) E Aip^ is a. minimizer for the problem ([7])-([8]), then the 
(y9j's solve (fTOj) with ei < 62 < ■ ■ ■ < Bn < the N lowest eigenvalues 
of the operator h^p [5]. 

In [5] it was proved that solutions {ipi, . . . ,(Pn} to — and, more 
generally, all eigenfunctions of the corresponding Hartree-Fock opera- 
tor h,p — are smooth away from x = (the singularity of V^), and that 
(for the (y9j's for which Ei < 0) they decay exponentially. (The solutions 
studied in [S] came from a minimizer of , but the proof trivially 
extends to the solutions {^J^jneN = {{¥^1 y ■ ■ ■ y V^n}} ® found 
in and to all the eigenfunctions of the corresponding Hartree-Fock 
operators mentioned above). The main theorem of this paper is the fol- 
lowing, which completely settles the question of regularity away from 
the origin of solutions to the equations ([9]). 

Theorem 1.1. Let Za < 2/tt, and let N >2 be a positive integer such 
that N < Z + 1. Let if = {(^1, . . . , (^tv}, e H^/^{R^), i = l,...,N, 
he solutions to the pseudorelativistic Hartree-Fock equations in ([9]). 
Then, for i = 1, . . . , N , 

if.eC^iR'MO}), (14) 

that is, the Hartree-Fock orbitals are real analytic away from the origin 
in R^. 
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Remark 1.2. (i) The restrictions < 2/7r, iV < Z + 1, and > 2 
are only made to ensure existence of i/^/^-solutions to ([9]). In fact, 
our proof proves analyticity away from x = for if^^^-solutions to (jH]) 
for any Za. For the case iV = 1, (|9]) reduces to {T — V)(p = eip and 
our result also holds for if^/^-solutions to this equation (see also (iv) 
and (v) below about more general V for which the result also holds for 
the linear equation). More interestingly, the result also holds for H^^^- 
solutions to ([1]) (which, strictly speaking, cannot be obtained from Q 
by any choice of A^). 

(ii) The statement also holds for any eigenfunction of the associated 
Hartree-Fock operator given by ( JTTi) . 

(iii) It is obvious from the proof that the theorem holds true if we 
include spin. 

(iv) As will also be clear from the proof, the statement of Theorem ll.il 
(appropriately modified) also holds for molecules. More explicitely, for 
a molecule with K nuclei of charges Zi, . . . , Zk, fixed at . . . , Rk € 

replace Vin^hj J2k=i with V^fc(x) = Zka/\x-Rk\, Zka < 2/%. 
Then, for A^ < 1 + Ylk=i Hartree-Fock minimizers exist (see 
Remark 1 (viii)]), and the corresponding Hartree-Fock orbitals are real 
analytic away from the positions of the nuclei, i.e., belong to C"^(M^ \ 
{Ri, . . . , Rk})- 

(v) Another approximation to the full quantum mechanical problem 
is the multiconfiguration self- consistent field method (MC-SCF). Here 
one minimizes the quadratic form q defined by the operator H given 
in ([3]) (or, more generally, with V from (iv)) over the set of finite sums 
of Slater determinants instead of only on single Slater determinants as 
in Hartree-Fock theory. If minimizers exist they satisfy what is called 
the multiconfiguration equations (MC equations). For more details, 
see [TOl [T3l 122] . As will be clear from the proof, the statement of 
Theorem 11.11 also holds for solutions to these equations. 

(vi) In fact, for V we only need the analyticity of V away from finitely 
many points in M'^, and certain integrability properties of Vipi in the 
vicinity of each of these points, and at infinity; for more details, see 
Remark 14.11 

(vii) As will be clear from the proof, the statement of Theorem 11.11 
also holds for other non-linearities than the Hartree-Fock term in ([H]), 
namely Iv^l^v^ as in ([2]) (for k even; for k odd, one needs to take ip^^^). 
The L^-space in which one needs to study the problem (see Proposi- 
tion 12.11 and the description of the proof below for details) needs to be 
chosen depending on k in this case (the larger the k, the larger the p). 
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(viii) Also, as will be clear from the proof, the result holds if T(— iV) = 
|V| (i.e., T(p) = IpI) in Q. In ^ below, ^(p)"^ should then be 
replaced by (|p| + 1)"^ (and '1' added to 'a~^ The only proper- 
ties of E{p)~^ used are in Lemmas IC.ll and IC.2| which follow also for 
(|p| + from the same methods with minor modifications. Similarly, 
one can replace T(p) with (—A + a~'^y, s G [1/2, 1]. 

(ix) The result of Theorem 11.11 in the non-relativistic case (T(— iV) 
replaced by —aA in ([3])) was proved in [131 EI]; see also the discussion 
below. In this case, it is furthermore known [TU] that, for x G Bj.{0) for 
some r > 0, (^^(x) = ipf\^) + \^\(pl^\^) with (pi^\(pf^ G C"(5,(0)). 

Combining the argument in [131 121] with the analyticity away from 
the position of the nucleus of solutions to the MC equations (see Re- 
mark [L2l (v)) we readily obtain the following result. 

Theorem 1.3. Let \E' 6e a (quantum) ground state of the operator 
H given in ([3]). Then ^ is not a finite linear combination of Slater 
determinants. 

Remark 1.4. The same holds with V as in Remark 11.21 (iv). 

Description of the proof of Theorem The proof of Theorem 11.11 is 
inspired by the standard Morrey-Nirenberg [2^ proof of analyticity of 
solutions to general (linear) elliptic partial differential equations with 
real analytic coefficients by 'nested balls'. A good presentation of this 
technique can be found in [16]. (Other proofs using a complexification 
of the coordinates also exist and have been applied to both linear and 
non-linear equations; see [26j and references therein.) 

In |16j one proves L^-bounds on derivatives of order k of the solution 
in a ball Br (of some radius r) around a given point. These bounds 
should behave suitably in k in order to make the Taylor series of the 
solution converge locally, thereby proving analyticity. 

The proof of these bounds is inductive. In fact, for some ball Bn with 
R > r, one proves the bounds on all balls Bp with r < p < R, with 
the appropriate (with respect to k) behaviour va. R — p. The induction 
basis is provided by standard elliptic estimates. In the induction step, 
one has to bound k + 1 derivatives of the solution in the ball Bp. To do 
so, one divides the difference Br \ Bp into k + 1 nested balls using k + 1 
localization functions with successively larger supports. Commuting 
m of the k derivatives (in the case of an operator of order m) with 
these localization functions produces (local) differential operators of 
order m — 1, with support in a larger ball. These local commutator 
terms are controlled by the induction hypothesis, since they contain 
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one derivative less. For the last term — the term where no commutators 
occur — one then uses the equation. 

This approach poses new technical difficulties in our case, due to the 
non-locality of the kinetic energy T(p) = ^/ — A + — and the 
non-linearity of the terms R^ipi and K^pipi. 

The non- locality of the operator V — A + a~'^ implies that, as op- 
posed to the case of a differential operator, the commutator of the 
kinetic energy with a localization function is not localized in the sup- 
port of the localization function. That is, when resorting to proving 
analyticity by differentiating the equation, the localization argument 
described above introduces commutators which are (non-local) pseudo- 
differential operators. Now the induction hypothesis does not provide 
control of these terms. Furthermore, it is far from obvious that the 
singularity of the potential V outside Br does not influence the regu- 
larity in Br of the solution through these operators (or rather, through 
the non- locality of V — A + Loosely speaking, the singularity of 

the nuclear potential 'can be felt everywhere'. (Note that if we would 
not have a (singular) potential V one could proceed as in [TT] and 
prove global analyticity by showing exponential decay of the solutions 
in Fourier space.) 

We overcome this problem by a new localization argument which 
enable us to capture in more detail the action of high order deriva- 
tives on nested balls (manifested in Lemma [B. II in Appendix [B] below) . 
This, together with very explicit bounds on the (smoothing) opera- 
tors (j)E{p)~^D'^x for X and with disjoint supports (see Lemma [Q2l) . 
are the main ingredients in solving the problem of nonlocality. The 
estimates are on (j)E{p)~^D'^x (iiot 'pE{p)D^x), since we invert -E(p) 
(turning the equation into an integral operator equation, see (l35l)). Our 
method of proof would also work in the non-relativistic case, since the 
integral operators (—A + 1)~^ and -E'(p)~^ enjoy similar properties. 

The second major obstacle is the (morally cubic) non-linearity of the 
terms Ripipi and K^ipi. 

To illustrate the problem, we discuss proving analyticity by the above 
method (local L^-estimates) for solutions u to the equation Au = u^. 
When differentiating this equation (and therefore m^), the application of 
Leibniz' rule introduces a sum of terms. After using Holder's inequality 
on each term (the product of three factors, each a number of derivatives 
on u), one needs to use a Sobolev inequality to 'get back down to L^' in 
order to use the induction hypothesis. Summing the many terms, the 
needed estimate does not come out (in fact, some Gevrey- regularity 
would follow, but not analyticity). 
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In the quadratic case this can be done (that is, for the equation 
Au = this problem does not occur), but in the cubic case, one 
looses too many derivatives. 

The second insight of our proof is that this problem of loss of deriva- 
tives may be overcome by characterizing analyticity by growth of deriva- 
tives in some with p > 2. When working in for p > 2, the loss of 
derivatives in the Sobolev inequality mentioned above is less (as seen 
in Theorem ID.ip . Choosing p sufficiently large allows us to prove the 
needed estimate. The operator estimates on (f)E{p)~^D^x mentioned 
above therefore have to be L^-estimates. In fact, using — L'^ esti- 
mates, one can also deal with the problem that the singularity of the 
nuclear potential V 'can be felt everywhere'. 

Note that taking p = oo would avoid using a Sobolev inequal- 
ity altogether (L°° being an algebra), but the needed estimates on 
(f)E{p)~^D'^X cannot hold in this case. For local equations an approach 
to handle the loss of derivatives (due to Sobolev inequalities) exists. 
This was carried out in [T2], where analyticity of solutions to elliptic 
partial differential equations with general analytic non-linearities was 
proved. Friedman works in spaces of continuous functions. In this ap- 
proach, one needs to have a sufficiently high degree of regularity of the 
solution beforehand (it is not proved along the way). Also, since the 
elliptic regularity in spaces of continuous functions have an inherent 
loss of derivative, one needs to work on a sufficiently small domain in 
order for the method to work. We prefer to work in Sobolev spaces 
since this is the natural setting for our equation and since the needed 
estimates on the resolvent are readily obtained in these spaces. 

For an alternative method of proof (one fixed localization function, to 
the power k, and estimating in a higher order Sobolev space (instead 
of in L^) which is also an algebra), see Kato [18] (for the equation 
Au = v^) and Hashimoto [TT| (for general second order non-linear 
analytic PDE's). 

Additional technical difficulties occur due to the fact that the cubic 
terms, -RtpV'i and K^ipi, are actually non-local. 

Note that in the proof that non-relativistic Hartree-Fock orbitals 
are analytic away from the positions of the nuclei (see [121 [22]), the 
non-linearities are dealt with by cleverly re-writing the Hartree-Fock 
equations as a system. One introduces new functions (pij = [^PiTp]] * 
I ■ |~^, which satisfy — A0jj = AiTifilpJ. This eliminates the terms 
R^ipi, K^LPi, turning these into quadratic products in the functions 
(fii, (pi J, hence one obtains a (quadratic and local) non- linear system of 
elliptic second order equations with coefficients analytic away from the 
positions of the nuclei. The result now follows from the results cited 
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above [IHIEE]. (In fact, this argument extends to solutions of the more 
general multiconfiguration self-consistent field equations, see [T31I22] .) 

This idea cannot readily be extended to our case. The operator E{p) 
is a pseudodifferential operator of first order, so when re-writing the 
Hartree-Fock equations as described above, one obtains a system of 
pseudodifferential equations. This system is, as before, of second (dif- 
ferential) order in the auxiliary functions (pij, but only of first (pseu- 
dodifferential) order in the original functions ipi. Hence, the leading 
(second) order matrix is singular elliptic. Hence (even if we ignore the 
fact that the square root is non-local) the above argument does not 
apply. 

To summarize, our approach is as follows. We invert the kinetic 
energy in the equation for the orbitals thereby obtaining an integral 
equation to which we apply successive differentiations. The localization 
argument of Lemma IB. II together with the smoothing estimates on 
(f)E{py^D(^X handle the non-locality of this equation. By working in 
LP for suitably large p one can afford the necessary loss of derivatives 
from using Sobolev inequalities when treating the non-linear terms. 

2. Proof of analyticity 

In order to prove that the (fi's are real analytic in \ {0} it is 
sufficient pO^ Proposition 2.2.10] to prove that for every xq G \ {0} 
there exists an open set f/ C \ {0} containing xq, and constants 
C,7^ > 0, such that 

\d^Vi{^)\ ^ ^ ^ for all X G f/ and all /3 G N;^ . (15) 

Let xo G M'^ \ {0}, and let u be the ball i?i?(xo) with center xq and 
radius R := min{l, |xo|/4}. For 5 > we denote by us the set of points 
in u at distance larger than 6 from du, i.e., 

:= {x G w I (i(x, du) > 5} . (16) 

By our choice of uj we have us = B^^si'^o)- Therefore cj^ = for 5 > R. 
In particular, by our choice of R, 

us = ^ for 5>l. (17) 

For Vt C M" and p > 1 we let L^(f2) denote the usual L^-space with 

norm = (/^ |/(x)|Pc?x)'/^ We write ||/||, ^ ||/||lp(m3). In 

the following we equip the Sobolev space iy"'^(f2), Q C M", m G N 
and p G [1, oo), with the norm 

||M||vK'"'P(n) := ^ WD^uWi^viji) . (18) 

|(T|<m 
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Theorem 11.11 follows from the following proposition. 

Proposition 2.1. Let Za < 2/tt, and let N > 2 be a positive integer 
such that N < Z + 1. Let cp = {ipi, . . .^^Pn}, e H'^^'^iM?), i = 
1, . . . ,N, be solutions to the pseudorelativistic Hartree-Fock equations 
in(Q. Lei xo G M^\{0}, i? = min{l, |xo|/4}, anc? w = i?i?,(xo). Define 
Us = Br_s{^o) for6>0. 

Then for all p > 5 there exist constants C,B > 1 such that for all 
j e N, for all e > such that ej < R/2, and for all i E {I, . . . , N} we 
have 

el^l||DVlUp(a;,,) < C^l'^l for all /3eNl with |/3|<j. (19) 

Given Proposition 12.11 the proof that the y^j's are real analytic is 
standard, using Sobolev embedding. We give the argument here for 
completeness. We then give the proof of Proposition 12.11 in the next 
section. 

Let U = i?R/2(xo) = coji/2 C CO. Using Theorem ID. 51 and (fT9l) we 
have (fi G C{U). Therefore it suffices to prove (JTSl) for \(3\ > 1. Fix 
i G {1,...,A^} and consider /3 G Nq \ {0} an arbitrary multiindex. 
Setting j = \I3\ and e = {R/2)/ j it follows from Proposition 12.11 (since 
ej = R/2) that there exists constants C, B > 1 such that 

l|flV.II„,„.„,<c(f)'^' = c(^)'^'|/J|l^ (20) 

with C, B independent of the choice of /3. By Theorem ID. 51 (see also 
Remark ID.6P there exists a constant = i^4(p, xq) such that, for all 
/3'eNl\ {0}, 

sup |i^^V.(x)| < ^4 V P^'+>.|Up(.,/,) 

1^1 

kl<i 

using ([20]). Using that R<1 < B, that #{(t G N|j | |a| = 1} = 3, and 
that, from flA.7p . 

(i + i/3i)^^"''<^e^|^'i mi. 

this implies that for all /3' G N|] \ {0}, 

sup|Z)^^.(x)|<(-^)(— ) l^'l!. (21) 
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Since |o-|! < Sl'^'a! for all a eNl (see flA.4p in Appendix S below), this 
implies that 

sup|D^V,,(x)| <C^, (22) 

for some C,TZ > 0. This proves f[T^ . Hence is real analytic in 
M'^ \ {0}. This finishes the proof of Theorem 11.11 
It therefore remains to prove Proposition 12.11 

Remark 2.2. We here give explicit choices for the constants C and B 
in Proposition 12.11 
Let 

\'Pa{y)My)\ 



Ci := max 

l<a,b<N 



■ -y| 



■dy 



(23) 



Note that by below, this is finite since G H^/'^{R^), i = l,...,N. 
Furthermore, let A = A(xo) > 1 be such that, for all a G Nq, 

sup|L'"y(x)| < yll"l+V|!. (24) 

xSo; 

The existence of A follows from the real analyticity in u = i?j?(xo) (re- 
call that R = min{l, |xo|/4}) of = Za\ ■ (see e. g. [20, Proposition 
2.2.10]). Assume without restriction that A > a'^ + maxi<j<jv 

Let Kx = Ki{p),K2 = K2{p), and = K^{p) be the constants 
in Lemma \C1\ Corollary ID.2[ and Corollary ID.4[ respectively (see 
Appendices [Cl and [D] below) . Then let 



C2 = max {7^1,256 V2/7r} , (25) 
C3 = max {47r(l + 2Ci/R^)K^, im-nKlK^] . (26) 

Choose 

C > max {1, ||v3i||iyi>pH, \Ml'^p{B2r{^o))^ —\^o\^^'^~^^^^'^^^M2, 

i£{l,...,N} vr 

.48^2, ^„ /-^ N 153672.,, „ . 

A + 48V2C1—+ ]yih}. (27) 

vr Ztc tt^IxqI 

That C < 00 follows from the smoothness away from x = of the yjj's 
Theorem 1 (ii)] and the fact that, since ipi e H^^'^(M:^), 1 <i <N, 
we have (fi G L^(]R^), 1 < i < A^, by Sobolev's inequality. Then choose 

1 6 

B > max {48AC2, a, — , ACl {lQ0C^K2C3f, {2ANC2lZ)\ IQK^] , 

|xo| 

(28) 
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where is the constant (related to a smooth partition of unity) intro- 
duced in f IB.Sp . In particular, B > 48. We will prove Proposition 12.11 
with these choices of C and B. 

3. Proof of the main estimate 

We first make ([6]) more precise, thereby also explaining the choice of 
A^TV in dH]). By Kato's inequality [HI (5.33) p. 307], 

/ i^dx<^ / \p\\f{p)\'dp for feH'/\R') (29) 

(where /(p) = (27r)^'^/^ J^^ e~'^'P/(x) dx denotes the Fourier transform 
of /), and the KLMN theorem [291 Theorem X.17] the operator ho given 
as 

ho = T(-iV) - V (30) 

is well-defined on H^^'^(R^) (and bounded below by — a^^) as a form 
sum when Za < 2/tt, that is, 

(m, V) = {E{py/\,E{py/\) - a-\u,v) - {V^/^u,V^/''v) 

for u,ve H^/^M.^) . (31) 

By abuse of notation, we write E{p) for the (strictly positive) operator 
£'(— iV) = V— A -|- q;~^. For {ipi, . . . ,ipiy) G A^at, the function 
given in f lT2|) belongs to L°°(M'^) (using Kato's inequality above), and 
the operator given in f|T3|) is Hilbert-Schmidt (see [5] Lemma 2]). 
As a consequence, when Za < 2/7r, the operator h^ in (fTTj) is a well- 
defined self-adjoint operator with quadratic form domain if^/^(M^) such 
that 

{u,h^v) = {u,hov)+a{u,R^v)-a{u,K^v) ioi u,v e H^^'^{R^) . (32) 

Since {u, Ripu) — (m, K^u) > for any u G L^(M^), also h^ is bounded 
from below by — a~^. 

Then, for (ui, . . . ,um) G AIat, the precise version of (j6]) becomes 

E^^{ui, ...,un) 

-- T [ [ ^^f^)^^.(y) (33) 

The considerations on i?,^ and i^',^ above imply that also the non-linear 
terms in ([33]) are finite for m G ^/^/^(M^), 1 < z < A. 
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If {ipi, . . . , y^Ar) G A^AT is a critical point of in ( l33l) . then (p = 
{ipi, . . . ,{p]y} satisfies tlie self-consistent HF-equations ffTU]) with the 
operator defined above. 

Note that E{p) is a bounded operator from H^^'^{R^) to H-'^^'^{R^), 
and recall that fl29|) shows that V also defines a bounded operator from 
i/i/2(M3) to if-V2(M3) (for any Za). As noted above, both i?<^ and 
are bounded operators on L^(R'^) when {ipi, . . . ,ipiy) E A4n- In 
particular, this shows that if {ipi, . . . , ip^) E J^n solves ffTOj) . then 

E(p)v9i - a" Vi - ^V^i + «-R<pV^i ~ ctK^ipi = Eiipi , I <i <N , (34) 

hold as equations in if~^/^(R'^). Using that £'(p)~^ is a bounded op- 
erator from i7~^/^(]R^) to if^/^(]R^), this implies that, as equalities in 
iJ^/^(]R^) (and therefore, in particular, in L^(]R^)), 

+ aE{p)-'K^ipi + {a-' +ei)E{p)-^ipi , l<i<N, (35) 

Proof of Proposition \2.1\ : The proof of Proposition 12. II is by induction 
on J G Nq. More precisely: 

Definition 3.1. For p > 1 and j G Nq, let V{p,j) be the statement: 
For all e > with ej < R/2, and alH G {1, . . . , A^} we have 

el^lp^</^^||L.K,) < C^l^l for all f3eNl with \f3\ < j , (36) 

with C, -B > 1 the constants in Remark I2.2[ 

Then Proposition 12.11 is equivalent to the statement: For all p > 5, 
V{p,j) holds for all j G Nq. This is the statement we will prove by 
induction on j G Nq. 

Induction start: For convenience, we prove the induction start for 
both j = and j = 1. 

Note that V{p,0) trivially holds since (see Remark 12.21) 

C = C{p) > max^ \\<Pi\\LP{uj) ■ (37) 

Also V{p, 1) holds by the choice of C, since 

C = C{p)> max \\D^ipi\\Lpr^) . (38) 

l<i<N, 

i/e{i,2,3} 

Namely, since C u, (1361) holds for = (and all e > 0) using (!37|) . 
For /3 G No with = 1 = j (i.e., f3 = for some u G {1,2,3}), and 
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all e > with e = ej < R/2 < 1, 

<C <CB = CB^^^ . (39) 

Here we again used that oue C uj, (l38i) . and that B > 1 (see Remark [2^ . 
We move on to the induction step. 

Induction hypothesis: 

Let p > 5 and j G No, j > 1. Then V{p,j) holds for all j < j. (40) 

We now prove that V{p,j + 1) holds. Note that to prove this, it 
suffices to study /3 G Nq with \f3\ = j + 1. Namely, assume e > is such 
that e(j + 1) < R/2 and let /3 G N|] with < j + 1. Then < j and 

< -R/2 so, by the definition of us and the induction hypothesis, 

el^l||/^^.l|L.K,,,,,) < 6l^l||D^.|U.K,) < CB\^\ . (41) 
It therefore remains to prove that 
el^l||DV,|Up(<.^^^^^^) < CB\^\ for all e > with e{j + 1) < R/2 

and all /3 G N|j with \/3\ = j + 1 . (42) 

Remark 3.2. To use the induction hypothesis in its entire strength, it 
is convenient to write, for £ > 0, e > such that ei < R/2, and a G Nq 
with < |cr| < j, 

ei - 

WD^'ipiWLPiui.e) = WD'^ipiWiPiuj--.) with e = j = \a\ , 

^■^ \(7 \ 

SO that, by the induction hypothesis (applied on the term with e and 
j) we get that 

< =c[LA^ (-) . (43) 

Compare this with (1361) . With the convention that 0° = 1, (l43i) also 
holds for |cr| = 0. 

We choose a function $ (depending on j) satisfying 

$ G C^iuJe{j+3/4)) , < $ < 1 , with $ = 1 on . (44) 

Then 

II^VllL.(a.,(,+„) < W'^D^VX- (45) 

The estimate (l42|) -and hence, by induction, the proof of Proposition l2.H — 
now follows from the equations (135|1 for the ipi^s, (l45l) and the following 
two lemmas. 
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Lemma 3.3. Assume ( l40l) (the induction hypothesis) holds. Let $ 
be as in flH|) . Then for all i G {1,...,A^}, all e > with e(j + 
1) < R/2, and all /3 G with \f3\ = j + I, both (!>Df^E{p)-^VyDi and 
$D^E(p)^Vi belong to Lp{R^), and 

\\^D^E{pr'V^^l<^[^f\ (46) 

||(«-^ + 5,,)$D''E(p)-V.||p < j(f , (47) 
where C,B > 1 are the constants in fl5Bl) (see also Remark \2.2^) . 

Lemma 3.4. Assume (HOl) (the induction hypothesis) holds. Let $ be 
as in (jBD. Then for all i e {1,. . .,N}, all e > with e(j + 1) < 
R/2, and all /3 G Njj with \f3\ = j + I, both <!>D^E{p)-^R^<^i and 
<^Dl^E{py^K^'^i belong to Lp{R^), and 

\\a^D^Eip)-'R^^,\\,<^(^^Y\ 

\\a^D^E{pr'K^^4,<^(^jY\ 

where C,B > 1 are the constants in ( |36l) (see also Remark \2.2\} . 
Remark 3.5. For a,b E {1, . . . , A^}, let Ua^b denote the function 



Ua,b{^) = / — . 7— ay, xgM . (48) 

In particular, ||f/a,6||oo < Ci for all a, 6 G {1, . . . ,N} (see (l23l) ). Note 
that (see and (fT^) 

i^cpV^i = ^ Ui^E^i , K^Lfi = ^ Ui^e^e . (49) 

£=1 £=1 

Hence Lemma 13.41 follows from the following lemma and the fact that 
Za < 2/71 < 1. 

Lemma 3.6. Assume fHOl) (the induction hypothesis) holds. Let $ be 
as in fj44|) . For a, 6 G {1, ... , A^}, /et Ua^ be given by fHHj) . T/ien /or 

a, 6, i G {1, . . . , A^}, all e > with e{j + 1) < i?/2, anc? all /3 G 
wt/i = j + 1, <l>Df^E{p)-^Uam belong to Lp{R^), and 

||$D^i?(p)-f/,,^.||p<^(f)'", (50) 
where C,B > 1 are the constants in (!36l) ('see also RemarklEM)- 
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It therefore remains to prove Lemmas 13.31 and 13.61 This will be done 
in the two following sections. □ 



4. Proof of Lemma [3731 
We prove Lemma [3.31 by proving (146!) and (1471) separately. 

Proof of din]) : Let (J e Mf] and u G {1, 2, 3} be such that P = a + e^, so 
that = DyD'^ . Notice that |cr| = j. Choose localization functions 
{Xfc}i=o a^'i {^fc}i=o a^ Appendix iBl below. Since V^pi G H~^/'^{M?), 
and ^(p)~^ maps //"(R^) to H'+^^R^) for all s G M, Lemma ED (with 
i = j) implies that 

j 

k=0 
k=0 

+ ^D,E{p)-'D'^[r],Vipi], (51) 

as an identity in H~^^^~^^^'^{M.^) (we have also used that -^(p)""^ com- 
mutes with derivatives on any H'{R^)). Here, [■ , ■] denotes the com- 
mutator. Also, |/3fc| = k, = 1, and < rjk, Xk < 1- (For the support 
properties of r]k,Xk, see the mentioned appendix.) We will prove that 
each term on the right side of (1511) belong to L'P(]R^), and bound their 
norms. The proof of (H^ will follow by summing these bounds. 

The first sum in (15T1) . Let 6k be the characteristic function of the sup- 
port of Xk (which is contained in u). Since V is smooth on the closure 
of u it follows from the induction hypothesis that the D"""^*: [Vipi] 's be- 
long to LP{u') for any u' CC u. Also, the operator ^Di,E{p)^^D'^''Xk 
is bounded on ^^(IR^) (as we will observe below). Therefore we can 
estimate, for A; G {0, . . . , j}, 

\\<l>D,E{pr'D''''XkD''-^''[Vip,]\\p 

= II imPr'D.D^'XkWkD'^-^'' [Vv^] lip 

< \\<^E{p)-'D,D^'xk\\B, \\ekD"-f'ny'P^]\\p- (52) 

Here, || ■ is the operator norm on Bp := B{Lp(M.^)), the bounded 
operators on Lp(R^). 
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For k = 0, the first factor on tlie riglit side of ( l52l) can be estimated 
using Lemma IC.l] (since \/3o\ = 0). This way, since ||xo||oo = ll'^lloo = 1, 

\\^E{p)-'D,xo\\b,<K,, (53) 

with Ki = Ki{p) the constant in (IC.ip . 

For A; > 0, the first factor on the right side of (1521) can be estimated 
using (IC.4I) in Lemma [Cl2] (with r = 1, q* = p = p). Since 

dist(supp Xki supp $) > t{k — 1 + 1/4) 

and llxfclloo = ll^lloo = 1, this gives (since (/3fc + ey)\ < {\(3k\ + 1)! = 
(fc + 1)!) that 

32V2{k + l] 



k Ve(A;- 1 + 1/4) 



<B^(?y. (54) 



IJBlpj-'D^C'Xtlls, < C2[-y , (55) 



TT V e. 

It follows from (ES]) and (El]) that, for all k e {0, . . . , j}, u e {1, 2, 3}, 

. e- 

with C2 as defined in (1251) . 

It remains to estimate the second factor in ( l52l) . Recall the definition 
of the constant A in dM]). It follows from ([MD and ([IT]) that, for all 
e > 0, £ e No, and cr G Ng, 

el"l sup |/^"l^(x)| < Al"l+V|! rl"l , (56) 

with C w as in defined in (fT6]) . 

For k = j, since /3j = cr, we find, by (156]) and the choice of C (see 
Remark that 

< CA . (57) 

The estimate for k G {0, ... ,j — 1} is a bit more involved. We get, 
by Leibniz's rule, that 

\\0kD''-^''[Vy:>,]\\p 

Now, supp 6k = supp Xk ^ uJe{j-k+i/A), so by ([56]), for all /i < a - 
ll^fc^'^^lloo < sup |D^V(x)| < e-I^UI''l+i|;u|!(j- A;)-!'^!. (59) 
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By the induction hypothesis (in the form discussed in Remark I3.2p . 
11/3, no"-/3fc-/i,^.|| < II r)°'^^'=^^(o. lU / 

^ ^ j^ k- /3fc-/i| j j k-&-Ml ^^^^ 

It follows from ( l58l) . ( I59|) . and ( l60l) that (using that |ct| = j, = k, 
and (lA.6p . summing over m = 

Note that, by (IATT]) . ioi <m < j - k, 

i — k\ m\{j — k — mV^^^"' e^^^'^Ji — k , , 

J \ \j ) ^ ^ ■> _ — < 1 _ (52) 



my (j — /i:)^ ^J j — k — m^^' 

To see the last inequality, look at the cases < m < (j — A;)/2 and 
j — k > m > {j — k)/2 separately. 

Hence (since B > 2A, see Remark 12. 2p . for any A; G {0, . . . , j — 1}, 

\\9,D--^^[V^,]\\, < CaQ'^'Y^^^Y < 2CaQ''\ (63) 

m=0 

Note that, by ( !57l) . the same estimate holds true if A; = j. 

So, from dig), (ESD, (El, the fact that e < 1 (since e{j + 1) < i?/2 < 
1/2), and the choice of B (in particular, B > 16; see Remark l2.2p . it 
follows that 



k=0 



^--^(f)'t(|)^^-(--^)(f)'-4(l)"' 

fc=0 

The second sum in ( 15T]) . Note first that [r]k,D'^''] = —{D'^'^qk) (recall 
that \fik\ = 1; see Lemma [B.ip . 

Comparing the second sum in fl^T]) with the first sum in flHTl) . one 
sees that the second sum is the first one with j replaced by j — 1 
and Xk replaced by —D^'^qk- Having now a derivative on the lo- 
calization functions we have one derivative less falling on the term 
V(fi. More precisely, the operator D''~^'^+^ contains \(j — /3a;+i| = 
j — {k + 1) = {j — 1) — k derivatives instead of |(7 — = j — k 
in D''~^'^. Then, to control D'^~^''+^\yLpi\ (with the same method 
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used above for D''~'^''[Vipi]) we need that supp D^'^r^fc is contained 
in a;e((j-i)-fc+i/4). Indeed we have much more: as for Xk we have 
suppL"''=?7fc C uj,^j_k+i/i) ^ We((i-i)-fe+i/4)- Finally, \\V'^'r]k\\oo < CJe, 
with > the constant in fIB.SP in Appendix [B] below. 

It follows that the second sum in (!5T|) can be estimated as the first 
one, up to one extra factor of C^/e and up to replacing j by j — 1 in 
the estimate (IMI) . Hence, using that e < 1, and the choice of B (see 
Remark I2.2p . we get that 

j-i 

'\^Ci,AC.^Q^-'<Ci,AC.)^^^^ (65) 

The last term in It remains to study 

<!>D^E{p)-'[r]^Vip,]. (66) 

We split V in two parts, one supported around x = 0, and one sup- 
ported away from x = 0, and study the two terms separately. We will 
prove below that this way, rjjVipi is actually a function in L^(M'^) + 
L^(M^). Upon using suitable operator bounds on ^D^E{p)~^x (for 
some suitable smooth x's)) combined with bounds on the norms of the 
two parts of rjjVipi, we will finish the proof. 

Let p = |xo|/4, and let 9p and 9p/2 be the characteristic functions of 
the balls -Bp(O) and -Bp/2(0), respectively. Choose Xp ^ C^(M^) with 
suppXp ^ -Bp(O), < Xp < 1, and Xp = 1 on 5^/2(0). Note that then 

dist(supp $, supp Xp) > = 2p , (67) 
by the choice oi u = Bu{yio), R = min{l, |xo|/4}, since supp $ C 

Now, 

^DPE{p)-%Vip,] = ^D^Eipr%Vxp^,] 

+ ^D^E{pr%V{l~XpWi]. (68) 

For the first term in ( |68|) . we use Lemma IC.2t with p = 1, q = 
p/{p — 1), and X = p. Then p, r G [1, oo) and q > 1, and q"^ + p^^ = 1. 
We get that (recall ( !67|l and that Xp^p = Xp): 

\\^D^Eip)-'[r]jVxpM\p < \\^D^E{pr%UJ\VjV9p^,h 



- ^^'(2^) i2pr/'-'{tm + 2)-5y\\Vep^^\\l■ (69) 
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Here we used that ||$||oo = ||Xp||oo = 1 and that r]j = 1 where 6p ^ 0. 
Note that j + I < e^^ (since, by assumption, e(j + 1) < R/2 < 1/2). 
Therefore, 

/3! < = (j + 1)! < (j + < e"^''+'^ = • (70) 
Note furthermore that since = j + 1 > 2 and t > 1, 

(t(|/3| + 2)-3)-^/^<l, (71) 
independently of (3. It follows that 
||4D''£(p)-'|%Vwlllp 

<i^(^)'"""||V.,,,|K(^)'r (72) 

Using Schwarz's inequality and that Za < 2/7r, 



ll^^'pV^illi < ll^^'plhll'/'ilb = ^aV|xo|7r||v5i||2 < ^ V |xo| Hv^iUs . (73) 



< 

TT 



(Note that WVOpUt < oo <^ t < 3.) It follows from ([72]), ([73]), and the 
choice of B and C (see Remark I2.2p that 

\\<S>D^E{p)-^lr,,Vx,'P,]\\p 

"^'x,|~.||,,|h(^)'"<|(f)'". (74) 

We now consider the second term in ( 168]) . Recall that $ is supported 
in C0e(j+i) and 

dist(supp $, supp r^j) > e(j + 1/4) . (75) 

Again, we use Lemma [C.2| this time with p = 3, q = p/{p — 1), and 
r = 3p/(2p+3). Then p^i + q-i+r"^ = 2, p G [1, cx)), q > 1, r G [1, 3/2) 
(since p > 3), and + p^^ = 1. This gives that 

||*r>''B(p)-'|>bV(l - XeMW, < \\<S>D'E{prh,,\\s,JV{l - x,)'P,h 

't^Kj(7^)""w^-^V4))''-(.(I/^|.2)-3)-'" 

X ll'^^(l -Xp)l|oo||¥'^||3• 
As before, we used that ||$||oo = ll^jlloo = 1- Note that 

/3!f^)'" < 321 ^ '^1' , = 32l^\ }' + ll < 32\^\ . (76) 
\j + 1/47 - (j + 1)1/^1 (j + 1)^+1 - ^ ' 



< 
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Since e(j + l) <R/2<1 and x < 3/2 it follows that (e(j + l/4))3A-2 < 
1. Also, by the choice of p, the definition of V , and since Za < 2/7r, 

|((l-^^p/2)l^)(x)| <^<^, xeM^ (77) 

|xo| 7r|xo| 

It follows from ((771) (and that < 1 - Xp < 1 -^p/2), (IH]), (176]), and the 
choice of C and S (see Remark 12.21) . that for alH = 1, . . . , iV (recall 
that = j + 1) 



4^2 16 „ „ /32\l/3| C 



<|.Z}^E(p)-i[r/,y(l-Xp)^.]||p 

^ II II 

TT 7r|Xo| V e / 24 V e / 
It follows from ([68]), ([71]), and (JTH]) that 

\\^dpe{^)-'[^,v^.^\1<^^[^Y\ ^^^^ 

The estimate (H^ now follows from (IFTj) and the estimates (IMj) . (ES]), 
and ([79]). □ 



Proof of Wl\ : Note that the constant functions Wj(x) = + Si 
trivially satisfies the conditions on V (= Za\ ■ \~^) needed in the 
proof above. In fact, having assumed A > + maxi<j<Ar \ei\ (See 
Remark [2?2D, ([24]) (and therefore ([56])) trivially holds for Wi. Also, 
for the term ^D'^ E{p)~^[7]jWiipi] we proceed directly as for the term 
E{p)~^[rijV {1 — Xp)fi] above (but without any splitting in Xp and 
1 - Xp), using that |PVi(x)| < A, x G M^. The proof of (jlTj) therefore 
follows from the proof of fl46p above, by the choice of C and B (see 
Remark 

This finishes the proof of Lemma 13.31 □ 

Remark 4.1. In fact, with a simple modification the arguments above 
(the local L^-bound on the two terms in (168]) ) can be made to work 
just assuming that, for all s > 0, 

V^, e L^BM) , V^, e L^R^ \ BM) ■ (80) 
5. Proof of Lemma [3T6] 

Proof of ( 150]) : Similarly to the case of the term with V in Lemma 13. 3[ 
we here use the localization functions introduced in Appendix IB] below. 
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With the notation as in the previous section (in particular, /3 = cr + 
with |cr| = j), Lemma [B. II (with i = j) imphes that 

j 

k=Q 

+ ^ <l>D,E{p)-'D''^-[r],, D^'^]D'^-''''+^[Ua,m] 

+ <!>D,E{p)-^D''[7]jUam]^ (81) 

as an identity in if~'^'(]R^). As in the proof of Lemma [■ , ■] de- 
notes the commutator, = k, \fik\ = 1, and < rjk, Xk < 1- (For the 
support properties of rik,Xk, see the mentioned appendix.) As in the 
previous section, we will prove that each term on the right side of flHTl) 
belong to L'f'(R^), and bound their norms. The claim of the lemma will 
follow by summing these bounds. 

The first sum in 0811) . We first proceed like for the similar sum in 
the proof of Lemma 13.31 (see flS^ . and after). Let 6^ be the charac- 
teristic function of the support of Xk- It follows from the induction 
hypothesis, using that —AUa,b = 4:HipaTpb, and Theorems ID. 51 and ID.3t 
that the ^'^"^''■[^a.feV'jl's belong to Lp{u') for any u' CC u. As be- 
fore, the operator ^Dt,E{j>)^^D^^Xk is bounded on ^^(IR^). Then, for 
/cG {0,...,j}, 

\\^D,E{p)-^D^>'XkD''-^''[Ua,m]\\p 

= II {^E{p)-^D,D^>'Xk)ekD''-^^ [UaM lip 

< \\^E{p)-^D,D^''Xkh^\0kD^-^'^[Ua,,^^]\\p. (82) 

The first factor on the right side of (!82|) was estimated in the proof 
of Lemma (see (ES])): For all A; e {0, . . . , j}, z/ e {1,2,3}, 

\\^E{p)-^D,DP^Xk\\B, < ^^ilY^ (83) 

with C2 the constant in (!25|) . 

It remains to estimate the second factor in (!82l) . For k = j, since Pj = 
a, we find that, by ( l23ll and the choice of C and B (see Remark 12. 2p . 

/5\l/2 ^ ^ 

\\OjUa,bVi\\p < \\Ua,b\\ 00 II V'i II LP(a;) 

<CiC<C(-) . (84) 
In the last inequality we also used that e < 1 (since e(j + l) < i?/2 < 1). 
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The estimate for /c G {0, . . . , j — 1} is more involved. We get, by 
Leibniz's rule, that 



< E ^')\\Ok{D^Ua,,){D'^'^''-^v^)\\p■ (85) 

We estimate separately each term on the right side of fl85|l . 
We separate into two cases. 

If /X = then, using the induction hypothesis (i.e., V{p,j — k); recall 
that supp6'fc C u^(^j_k)) and 

\\euUa,,D'^-P^^,\\,<C^ci^-) <l^[-) ■ (86) 

In the last inequality we used the choice of B (see Remark 12. 2p and 
that e < 1. 

If < /i < 0" — then (since suppxfc ^ ^aij-k+i/A)) Holder's 
inequality (with 1/p = l/(3p) + 2/(3p)) and Corollary ID. 21 give that 

\\ek{D''Ua,b){D''-^'-''Vi)\\p 

< \\ekD^Ua,bhp/2 \\0kD^-^''''V^hp 

II V"«lllV^'P(aj,(j_fc+i/4))ll V J 

Here, K2 is the constant in Corollary ID.2t and 6 = 2/p < 1. Note 
that a;e(j-A;+i/4) = -Br(xo) with r G [-R/2, 1], since e(j + 1) < R/2 and 
= min{l, |xo|/4} 

We will use Lemma [5.31 below to bound the first factor in flHTj) . The 
last two factors we now bound using the induction hypothesis. 

If G Nq is such that < ^ < a — (3^, then the induction hypothesis 
(in the form discussed in Remark 13. 2p gives (recall here f[T5]) and that 
W\ = jAf^k\ = k) that for the last two factors in (187|) we have 

||£)cr-/3fc-/i ||1-6I 

II ^*lli*'K(j-fc + l/4)) 



< 



'fj — k— \ i-fc-lMl /5\ i-fc-lMl] 1-6 
, \j -k + 1/4) [7) 
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and (using that B > 1 (see Remark [22!) e(j-^ + l/4) < e(j + l) < 
R/2 < 1) 



('^e(j-fc + l/4)) 



f j-k-\fi\ y-k-M/By- 
\j~k + l/4) \e) 



+ 3C 



j -k-\fi\ + ly-fc-lMl+i /5y-fc-lMl+i 
J -A; + 1/4 ) \~) 



< 



AC 



i - k-\ij,\ + ly-fc-lMl+Y^y-fc-lA^I+i 



j -k + 1/4 

It follows from (EH]) and ([89]) that for all /i G N[] with Q < ^ < a - 



9) 



'^«llLP(a;,(,_fc+i/4)) 



< C4'^(^) 



J - A; + 1/4 



(90) 



From (I57)) . Lemma [375] and (MI]) (using (1A.6P in Appendix \K\ below, 
summing over m = |/i|), it follows that 



0</i<CT-/3fe ^ 

m=l 



||^^,(OT„,,)(D'^-^'=-'^<^,)||p < c'ckJ-) X 



V e / 



j — A;\ (j — A; — m + 1 ) 



j—k—m+6 



> + 1/4)' 



X 



(j - A; + l/4)i-*^+« 
/ 1 , — / i?(rr 



■X 



/ 5(m+ 1/4) 
(j- A: + 1/4) J 



(91) 



Here, C3 is the constant from (I26|) . Recall also that 9 = 2/p. 
We prove that for m G {1, . . . , j — A;}, 



j — A;\ (j — A; — m + 1) 



j~k—m+9 1 



m + 1/4)' 



(j - A; + l/4)J^-'=+« 



< lOe-i/2+e 



Note first that, since e(j - A; + 1/4) < e(j + 1) < 1, 
(^■_fc + 1/4)1/2-^ <e-l/2+^ 



(92) 



(93) 



This shows that the inequality in ([92]) is true for m = j — k > 0, since 
< 1. For m < j — k, we use (IA.8I) in Appendix [A] below, and ( l93]) . to 
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get that (since (1 + 1/n)" < e) 

-k\{j -k-m + l)i-'=-'^+«(m + 1/4)'" 



m J {j -k + 1 /Ay-''+^ 

e25A2 _ A; - m + 1)^ I 
^/27^ [j-k-myl^ ^/m 

Since 9 <l/2 and m < j — A; — 1, we have that 



{j — k — m)i/2 



(95) 



The estimate f l92p for m G — A; — 1} now follows from 

(since 4^5/12/7? < 10). 

Inserting ( 192|) in ( 19T|) (and using again e(j — /c + 1/4) < 1 and 26' — 2 < 
0) we find that 



0<n<(7-l3k 



.ioc3C3.,(f)-n— i:[(j,)%^-^ 

m=l V ^ 



where we used that 9 < 2/5, B > A (see Remark 12.21) . and Ylm=i ^ 
1 + J-j°° (ix = 6 to estimate 

^/^X"^ 2 1 ^ 1 6 

m=l ' ' m=l * 

This is the very essential reason for needing p > 5. 
By the choice of B (see Remark 12.21) it follows that 

E (" : ll''.(fl'-C..)(fl°-*>-V.)ll,. < f(f )'"""' . (98) 

From (IH5]), (EE]), and (EH]) it follows that for all e {0, . . . , j - 1}, 

||6^,D'^-^'=[t/,,,(/.,]||p<C(-j . (99) 
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Using (IH2]), dHl, (El, and (EHD it follows for the first sum in (EI 
that 



J2 <^D,E{p)-'D^''XkD''-^' [Ua. 



k=0 



< 



fc=0 



k=0 



Since B > 16 (see Remark l2.2p the last sum is less than 2 and so for the 
first term in flHT]) we finally get, by the choice of B (see Remark 12. 2p 
that 



J2 ^B,E{p)-^D^'^XkB''-^'^[Ua,m 



k=0 



< ICoC 



B\j+y'^ CZ 



< 



12A^ V e / 



(101) 



The second sum in fl8T|) . By the same arguments as for the second sum 
in fl5T]) (see after fl64|) ). it follows that the second sum in f lHTj) can be 
estimated as the first one, up to one extra factor of C^/e (with > 
the constant in (IB.3P in Appendix |B] below) and up to replacing j by 
j — 1 in the estimate fllOip . Hence, by the choice of B (see Remark l2.2p 



J2'^B,E{p)-'D^''[r]k,D^'']D''-^''+'[Ua., 



k=0 



- e 12N \e) " 12A^ V e / 



(102) 



The last term in (ISTl) . Since a + Cu = (3, the last term in (ISTll equals 

^D^E{p)-'[r],Uam]- 

We proceed exactly as for the term ^D'^E{p)~^[rijV{l — Xp)fi] iii (pSj) 
(but without any splitting in Xp and 1 — Xp), except that the estimate 
in (jTfl) is replaced by ||f/a,6||oo < Ci (see (j23l) ). It follows, from the 
choice of B and C (see Remark 12. 2p that (recall that = j + 1) 

\\^D^E{p)-'[7],Ua,m]\\p < II'^'^''^(P)"'^.II% J|f/a,6</^.||3 



<r II 

^ <--l||¥'i||3| 



TT 



, , ■ (103) 

The estimate (150|) now follows from (IHTjl and the estimates ( llOip . fll02p . 
and ffT03|) . 
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This finishes the proof of Lemma 13.61 □ 



It remains to prove Lemma 15.31 below (L^^/^-bound on derivatives of 
the Newton potential Ua,b of products of orbitals, (pa^Pb)- 

In the next lemma we first give an L^^^^-estimate on the derivatives 
of the product of the orbitals (fi, needed for the proof of the bound in 
Lemma 15.31 below. 

Lemma 5.1. Assume ( HOl) (the induction hypothesis) holds. Then, for 
all a,b e {1, . . . , A^}, all /3 G N|] with < j - 1, and all e > with 
em + l)<R/2, 

, / R\ 1^1+261 

\\D^{VamL^./^i.^,,,,,,) < lOKlC^l + VW\)[-) , (104) 

with K2 from Corollaru W.^ C from Remark \2.S[ and 9 = 9[p) = 2/p. 
Proof. By Leibniz's rule and the Cauchy-Schwarz inequality we get that 

We use Corollary ID. 21 (with C0e{\i3\+i) = -Br(xo), r = R — e{\f3\ + 1); note 
that r e [R/2, 1], since e(|/3| + 1) < R/2 and R = min{l, |xo|/4}). This 
gives that, with K2 from Corollary ID. 2[ and 6 = 2/p, 

We now use the induction hypothesis (in the form discussed in Re- 
mark 13. 2p on each of the four factors in the sum on the right side 
of (fT05|) . Note that, by assumption, e(|/3| + 1) < ej < R/2 and 
< + 1 < + 1 < j (similarly, |/3 - /i| < |/3 - /i| + 1 < j). 
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Recalling ( |T8|) . we therefore get that, for all /z G Nq such that /i. < /3, 



< 



c 



- 



X 



C 



1/31 + 17 \e) 



1/31 + 1 

) 



- +3C 



e / 



1 



since (recall that e(|/3| + 1) < i?/2 < 1 and > 1) 



(|;,| + 1)1.1+1 -'- + 1)^" 

Proceeding similarly for the other two factors in (llOSp . we get (using 
f lA.6p in Appendix R] and summing over m = that 



5Z ( jll^VJU3PK(|^|+,))||/?'' ''<^b|U^''K(i;3i+i)) 



/^N 1/31+29 



<16''(CK2)2(^-j 

^ /|/3|\ [(m + ir+^ 

m=0 ^ ^ 



X 



(106) 



+ 1)1^1-™+!]^ [m^ 



^)|/3|-m]l- 



+ 1)1/31+29 



We simplify the sum in m. Note that for m = and m = the 
summand is bounded by 1. Therefore, for < 1 the estimate f ll04p 
follows from (1106^ . since 2 ■ 16^ < 7. It remains to consider > 2. 
For m > 1, m < we can use ( lA.Sp in Appendix lAl to get (since 
(1 + 1/n)" < e) that 



''V5a||L3p(a;,(|^l+l)) i'^*'('^6(|/3| + l)) 

A/27r V e / 



0<At</3 



X 



(1/31 + 1)1/31+29 

[(m + l)(|/3|-m + l)]' 



m=l 



Since the function 



/(x) = (x + l)(|/3|-x + l), xG [1, 1/31-1], 
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has its maximum (which is (|/3|/2 + 1)^) at x = |/3|/2, and since 

V 1 rm 1 

2^ ^ = < / , dx = TT, 



we get that 

0<M</3 ^'^^ 

< eV-(16e^)^y|(OT,)^v^(|)'"^''. (107) 

The estimate fll04p now follows from (11051) . (I106p . and (11071) . since (as 
P>5), 

eVi2(i6e2)ey|<^Q^ 2-16^<7. 

This finishes the proof of Lemma 15.11 □ 

The next two lemmas, used in the proof above of Lemma ESI control 
the L^P/^-norm of derivatives of Ua,b- 

Lemma 5.2. Define Ua.b by (148|) . Then for all a,b & {1, . . . , A^}, and 
all /i e N|j with < 2, 

WDf^UaAL^p/H.) < ^7rKs{C^ + 2Ci/R^) , (108) 

with K3 from Corollary D.4 , C from Remark \2.Sl Ci from ( l23l) . and 
R = min{l, |xo|/4}. 

Proof: Recall that u = Bpti^xo), R = min{l, |xo|/4}. Using (fTSj) . and 
Corollary ID. 4t we get that, for all /i G Nq with \fi\ < 2, 

\\D^Ua,b\\L3p/2{uj) < l|f^a,b||vy2,3p/2(B^(x(,)) (109) 

< Ks{\\AUa,b\\L--^P/2{B2R{^o)) + ;^ll^'i.felli3p/2(i32fl(xo))} • 

By the definition of Ua^b (see fHS]) ) we have 

- Af/,,fe(x) =47r(/^,(x)^(x) for x G M^ (110) 

and ||f/a,fe||oo < Ci (see (!23|) ). Hence, from fllOQp . Holder's inequality, 
and the choice of C (see Remark 12. 2( recall also that p > 5) 

||-D^t/a,b||L3p/2(^) < 47rK3{||v9a||i3p(B2^(xo))||v5fe||L3p(B2fl(xo)) 

+ -^||f/a,4oo|52/?(xo)r/^^} 

< 4717^3(^2 + 2Ci/R^) . 
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This finishes the proof of the lemma. □ 

Lemma 5.3. Assume fHOj) (the induction hypothesis) holds, and define 
Ua,b by m- 

Then for all a, 6 G {1, ... , A^}, all k e {0, . . . J - I}, all /i G N|] with 
l/^l ^ 3 — k, and all e > with e{j + 1) < R/2, 

\\D f/a,.|U3./.(.,,_,,,/,,) < CsC [j-k + l/J 

.c3.w(f)'^'^^^-^(-^)'^'^^^-^ (111) 

with 9 = 9{p) = 2/p, C and B from Remark IK^ and C3 the constant 



m 

Proof: If m := < 2, fillip follows from Lemma 15.21 and the defi- 
nition of C3 in ([26]), since e(j - A; + 1/4) < e(j + 1) < R/2 < 1, and 
C,B>1 (see Remark 

If m := > 3 then we write fi = fim~2 + + e^^ with z/j G 
{l,2,3},i = 1,2, |/im-2| = m — 2. Then by the definition of the 
ly^'^P/^-norm (recall f|T8l) ) we find that 

with ei such that 

ei(m-l + l/4) = e(j-A; + l/4). (113) 

To estimate the norm in flll2p we will again use that f/^ 5 satisfies 
(111 op . Applying /}'^™-2 to (IllOp and using the elliptic a priori estimate 
in Corollary ID. 41 (with r = ri = i? — ei(m — 1 + 1/4) and 5 = 5i = ei/4; 
recall that Up = i?ij_p(xo)) we get that 

P"f/a,.|L3./2K(^_,^,/,,) < 47rir3||/^^--HV^aW)llL3./2(.,^(_,,) 

+ f/,.|L3./2(.^,^^_^,), (114) 



with = K^Ip) the constant in (ID.9p . Notice that for this estimate 



we needed to enlarge the domain, taking the ball with a radius ei/4 
larger. 

We now iterate the procedure (on the second term on the right side 
of (dUD), with Q (i = 2, . . . , [f J) such that 

iiim - 2z + 1 + 1/4) = e,_i(m - 2{t - 1) + 1) , (115) 
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and with r = Vi = R — ei{m — 2i + 1 + 1 / 4) and 6 = 6i = ii/A. Note 
that ffTT^ and ffTT^ imply that, for i = 2, . . . , [f J , 



> ej„i > . . . > ei = e 



j -k + 1/4 
m - 1 + 1/4 



and 



(116) 



ei(m - 2ii + 1) < ii-i{m - 2{i - 1) + 1) 

< . . . < h{m - 1) < e(j - A; + 1/4) 



(117) 



We get that (with n£=i = 1 |/im-2j| = m — 2i), 



< 47rir3X] [p'^-^^(¥'aW)llL3p/2(., 



j-1 



1=1 



T\(— 

'ii{m-2i + l)) 11 \ g2 



LfJ 



+ n 



16^ 

el 



|I^'^-^-^f/.,IL3./.(.^.^^^^_^^^^^^). (118) 



Using f lll6p . and Lemma 15.11 for each i = 1, . . . , [yj fixed (note that 
iiim - 2i + 1) < R/2 by ffTTTD since e(j + 1) < R/2) we get that 



^=1 



16/^3 



<2mlc^^{^^)' 



B\ m+2e-2^m - 1 + 1/4x^+26-2/16/1-3^-1 



with K2 from Corollary ID. 21 and 9 = ^(p) = 2/p. Here we also used 
that 1 + Vm-2i < 2^. Note that J2i=i {^QK^/ B"^)'-^ < 2 since 
> 32K3 (see Remark O). It foUows that 



i-l 



1=1 



i=i 



<m.K?.K,C^V^{^) (-^^) . (120) 



We now estimate the last term in flllSp . Let 6 = m — 2[yJ e {0, 1} 
(depending on whether m is even or odd). Then, using (1116^ and 
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Lemma [521 we get that 

I — I 



[n 



1 



a,fe|lL3p/2K-|^l(„_2L™j+i)) 



^ / e(j-fc + l/4) y 
V m - 1 + 1/4 / 

.4.A„:.2C,.^)c^(f)-(-^)". (m) 

Here we also used that m > 3 and f^s > 1 (See Corollary ID .41) . that 
C > 1 and 5 > IQK^ (see Remark and that e(j - A; + 1/4) < 1. 

Combining f lTTSl) . ffM . and ffT2T]) finishes the proof of ffTTTj) in the 
case m = |yu| > 3. 

This finishes the proof of Lemma 15.31 □ 



Appendix A. Multiindices and Stirling's Formula 

We denote Nq = N U {0}. For a = (cti, 0-2, 0-3) e N|j we let \a\ : = 
+ (72 + (J3, and 

D'' ■.= Dl'Dl'Dl\ D,:= -i^=: -id, , z/ = l,2,3. (A.l) 
This way, 

We let a! := ai!o'2!cr3!, and, for n G No, 



(A.2) 



With this notation we have the multinomial formula, for x = (xi, X2, X3) G 
and n G Nq, 

(a;i+a;2 + X3r = x'^. (A.3) 

Here, x'^ := x'^^x'^^x'^^ . It follows that 

|a|! < 3l''la! for all a e^l, (A.4) 
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since, using ( lAl2ll . that (1, 1, 1)^ = 1 for all n G N|], and ([Qjl . 



We also define 

for cr, /i G Ng with fj. < a, that is, < a^, v = 1,2,3. Note that for 
all cr e and /c G No (see [HI Proposition 2.1]), 

Finally, by [1, 6.1.38], we have the following generalization of Stirling's 
Formula: For m G N, 

1 

ml = v^27rm'"^2 exp(— mH ) for some = -(^fm) G (0, 1) , (A. 7) 

12m 

and so for n,m & N, m < n, 

13(71) d{m) 'd{n — m) 



m) 727rm™+i/2(ri - m)"-'"+V2 ^^P*^ 12r2 12m 12(n - m 



gi/12 



< ( \ i 



/27r m'"+i/2(n - 
Appendix B. Choice of the localization 

Recall that, for xq G \ {0} and R = min{l, |xo|/4}, we have 
defined u = i?i?(xo), us = -B_r-5(xo), and that e > is such that 
e{j + 1) < R/2. Also, recall (see dH])) that we have chosen a function 
$ (depending on j) satisfying 

$ G C^iuJe{j+3/i)) , < $ < 1 , with $ = 1 on WeO+l) . (B.l) 

For J G N we choose functions {Xk}i=o, and {r^^j-^^Q (all depending 
on j) with the following properties (for an illustration, see figures 1 and 
2). The functions {Xk}i=o are such that 

Xo e C^(w,(j+i/4)) with xo = 1 on uje{j+i/2) , 
and, for A; = 1, . . . , j, 

Xk GC^(We(j-fe+l/4)) 

with <; = ^ We(i-fc+l/2) \ C^e(j-fe+l + l/4) , 

Xk = on \ (We(j-fc+l/4) \ l^e(j-fe+l + l/2)) , 



34 A. DALL'ACQUA, S. FOURNAIS, T. 0. S0RENSEN, AND E. STOCKMEYER 
Finally, the functions {r]k}l.^Q are such that for k = 0, . . . , j , 

rik e C-(R^) with {"^'^l on u^^-k+i/A) , 

Moreover we ask that 
Xo + ?7o = 1 on 

Xfc + % = 1 on y a;,Q_fc+i+i/4) for k = 1, . . . J , (B.2) 

Vk = Xk+i + Vk+i on for /c = 0, . . . , j - 1 . 

Furthermore, we choose these localization functions such that, for a 
constant > (independent of e,k,j,(3) and for all /3 e Nq with 
\/3\ = 1, we have that 



\D^Xk{^)\ < 


^ and |DV(x)|<^, 


(B.3) 


for A; = 0, ... , j, and all x e 






oj = Bh(xo) 






1 


1 1 1 1 

R/2 e{j + l) ■■■ 2e e 





Figure 1. The geometry of a; = Bu{^o) and the Uek — -Bfl_efc(xo)'s. 




Figure 2. The locahzation functions. 
The next lemma shows how to use these localization functions. 



ANALYTICITY OF PSEUDORELATIVISTIC HARTREE-FOCK ORBITALS 35 

Lemma B.l. For j G N fixed, choose functions {Xk}i=o, and {rik}l^Q 
as above, and let cr e Ng with \a\ = j . For £ e N with i < j , choose 
multiindices {l3kYk=o ^^^^ that: 

\/3k\ = k for k = 0, . . . ,i , (5k-i < h for k = l,...,i, and /3i < a . 

Then for all g G S'{R^), 

e 

D'^g = Y,D^'XkD''~P'^g (B.4) 

A:=0 
£-1 

fc=0 

with fik = - /3k for k = 0, . . . ,i - 1 (hence, = 1). 

Proof. We prove the lemma by induction on i from i = 1 to i = j. We 
start by proving the claim for i = 1. By using property (1B.2|) of the 
localization functions and that /3i = /3o + /^o = /^o (since f3o = 0) we 
find that 

D'^g = xoD'^g + r^oD'^g = xoD^g + VoD^-^'+^'g . (B.5) 

The first term on the right side of fIB.Sp is the term corresponding to 
= in the first sum in flB.4p . In the second term in fIB.Sp . commuting 
the derivative through ?7o, we find that 

^^^^~/3i+Mo^ = D'''r]oD^-^^g + [7]o, D^'']D^-^' g . 
Since rjo = xi + Vi by property flB.2p . this implies that 

r]oD''-'^'+^'»g 

= D^^XiD^^^^g + D^^'qiD^^^^g + [r/o, D^"]D''-^^g. (B.6) 

The identity flRi]) for £ = 1 follows from (|R5|) and flRel) . 

We now assume that (lB.4p holds for £ — 1 for some £ > 2, i.e., 

£-1 

D'^g = Y,D^'XkD''-P^g (B.7) 

fc=0 
£-2 

k=0 

and prove it then holds for i. Since = (3i — fii-i we can rewrite 
the last term on the right side of flB.7p as 
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Again, commuting the /i^_i-derivative through rji^i this imphes that 

= D^'iT], + Xi)D''-^'g + D^^-^ [r7,_i, D^'-^]D''-^'g , (B.8) 

using (!R2l) . Collecting together (!R7l) and (iRS!) proves that (iRill holds 
for £. 

The claim of the lemma then follows by induction. □ 

Appendix C. Norms of some operators on Lp(R^) 
In this section we prove two lemmas on bounds on certain operators 



a 



-2 



involving the operator E{p) = \/—A + 

Lemma C.l. Let the operators Sy = E{p)^^Dy, v G {1,2,3}, he 
defined for f e S{R^) by 

(^./)(x) = (27r)-3/2 / ^-■^E{p)-'pJ{p)dp, 
with /(p) = (27r)~'^/^ Jjg3 e~'^'P/(x) the Fourier transform of f . (Here, 

P = {PUP2,P3)-) 

Then, for all p G (l,oo), S^, extend to bounded operators, S^, : 
LP(R3) ^ LP(M3), e {1,2,3}. Clearly, = \\S,\\t3„ 7^ f^- 

We let 

K,=Ki{p) := \\Si\\b,. (C.l) 

Proof. This follows from [301 Theorem 0.2.6] and the Remarks right 
after it. In fact, since (by induction), 

Dl{p,E{pr') = P,,.(p)E(p)-i-2H , 7 G Nji , 

for some polynomials P^,u of degree I7I + 1, the functions m^(p) = 
Pi,E{p)~^ are smooth and satisfy the estimates 

|D^m,(p)| <C^,,|p|-l^l, jeNl 

for some constants C^^i, > 0, which is what is needed in the reference 
above. □ 

For p, q G [1, 00], denote by || ■ H^p^ the operator norm on bounded 
operators from LP(R3) to L'^{R^). 

Lemma C.2. Forallp,x G [I, 00), q G (l,oo), withp^^ + q^^+x^^ = 2, 
all a>0, all f3 e N|] (with \f3\ > 1 if x = 1), and all <l>,xe C^iR^) n 
L°°(M3) with 

dist(supp(x), supp($)) > d , (C.2) 
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the operator ^E{p)-^Df^x is bounded from (R^) to (/."(R^))' = L''*(]R3) 
(with q^^ + q*^^ = Ij, and 

||4B(p)-'BMle„. (C.3) 

< ^^/?!Q'V-(r(|/!| + 2) -3)-"|H.|U||xlU. 

In particular, (when t = 1, i.e., q* = p), 



TT \l3\-l\d 
for all 13 e nl with > 1. 

Proof. We use duality. Let f,ge Note that, since $/, D^{xg) G 

L^(M^), the spectral theorem, and the formula 

—=, X > 0, (C.5) 



a/x TT Jq X^t ^t 

imply that 



1 df 

if, ^Eipr'D^xg) = - -^if, $ (-A + a-' + t)-'D^x9) ■ 

Jo Vt 

By using the formula for the kernel of the operator (—A + a^^ + 1)^^ 
(IX.30)], and integrating by parts, we get that 



{f,<i>E{pr'D^X9) 



= -/ 4/ /W'^W / -[D^{x9)]{y)d^dy 

Jo VtJu^ Jr^ 47r|x-y| 

dt r r / p-Va-2+t|x-yk 

= ^^/ 4/ /(x)$(x)/ (dJ^— -)x(yMy)rfxrfy. 

Jo VtV 47r|x-y| / 

Notice that the integrand is different from zero only for |x — y| > d, 
due to the assumption (1C.2I) . Hence, by Fubini's theorem, 

{f,<!?E{p)-'D^X9)= [ / F(x)i7(x-y)G(y)rfxdy, (C.6) 



with F(x) = /(x)<l>(x), G(y) = x{y)9{y), and 
H{z) = H^,^,d{z) 
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Now, by ( 1C.9I) in Lemma [C.3I below, uniformly for a > 0, 



ii^(z)i<i{|.|>,}(z: 

= M\-\>d}{z, 



An"^ |z| V |z| 
V2 /3\ / 8 



-v^|z|/2 ^ 

Vi 



TT^ |z| 



and so, for all a > 0, r G [1, oo), and all /3 G Nq (with > 1 if r = 1), 

lA 



From this, (]C.6p . and Young's inequality [211 Theorem 4.2] (notice that 
Cy < 1), follows that, with p, q, r G [1, oo), p"^ + q"^ + = 2, 

\{f,mPr'D^X9)\<\\FUHUGl 

<(4vr)V^^/3!Q'"ci3A-2(,(|^|+2)-3)-^/^||FlU|G||, 



< 



TT 



-/3!(J)'^'rf'/^-'(t(|/3| + 2)-3)-^/^||$|| 



oo II A II C!0 



P • 



Since iS(M^) is dense in both (R^) and L'^ (M^), this finishes the proof 
of the lemma. □ 



Lemma C.3. For all s > 0, ^eR^\ {0}, and (3 G Njj, 

V2f3\f 8 



OP — 

X I ,1 



< 



s|x| 



< 



/ 8 



„-s|x|/2 



(C.8) 
(C.9) 



Proof. We will use the Cauchy inequalities [TTl Theorem 2.2.7]. To 
avoid confusion with the Euclidean norm | ■ | (in or in C^), we 
denote by | ■ |c the absolut value in C 
Let, for w = {wi, W2, w^) G and r > 0, 

p3(w) = {z e \ \z, - w,\c < r , u = 1,2, 3} (C.IO) 

be the poly-disc with poly-radius r = (r, r, r). The Cauchy inequalities 
then state that if u is analytic in -P^ (w) and if sup2gp3(w) |'w(z)|c < M, 
then 



\d^u{w)\c < Ml3\ r-l^l for all l3eNl. 



(C.ll) 
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We take w = x G M'^ \ {0} C and choose r = |x|/8. We prove below 
that then we have (with t? := ^21=1 ^ ^) 

Re(z2) > i|x|2 for z G P,^(x) . (C.12) 

It follows that a/z^ := exp(|Logz^) is well-defined and analytic on 
Pr{^) with Log being the principal branch of the logarithm. 
We will also argue below that 

1 



Re(Vz2) > ^|x| for z G P^(x) . (C.13) 
Then (by flal2D ) for all z G P^ip^), 

Iv^lc = > > |x|/V2 , (C.14) 

and (by (laTSj) ). for all s > and all z G P^^(x), 

I exp(-sv^)|c = exp(-sRe(v^)) < exp(-s|x|/2) . (C.15) 

Therefore, (^M) and dClQ]) fo llow from dCH]), (IClil) . and (lOTsD . 
It remains to prove (1C.12P and (IC.lSp . 

For z G P^''(x), write z = x + a + ib with a, b G satisfying 
\z,-x,\l = al + bl<{\^\/8y. Then 

z^ = |x + ap - |bp + 2i(x + a) ■ b , 

so, with e = 1/8, 

Re(z^) = |xp + |ap + 2x - a - |bp 

> (l-e)|x|2 + (2-e-i)|ap-(|a|2 + |b|2) 

> — xp > - xP . 

- 64' ' 2' ' 

This establishes (IC.12P . 

It follows from flC.12p that, with Arg the principal branch of the 
argument, 

- ^ < lArg(z2) < ^ for z G P^(x) . (C.16) 

Furthermore (still for z G P^^(x)), because of ( 1C.16I) . 

Re(v^) = |z2|^/'cos(iArg(z2)) > \z^\l/^/V2. (C.17) 

Combining with flcn2l) we get flClTal) . 

This finishes the proof of the lemma. □ 
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Appendix D. Needed results 

In this section we gather some results from the hterature which are 
needed in our proofs. 

Theorem D.l. [2, Theorem 5.8] Let Q be a domain in R" satisfying 
the cone condition. Let m G N, p G (1, oo). If mp > n, let p < q < oo; 
if mp = n, let p < q < oo; if mp < n, let p < q < p* = np/{n — 
mp). Then there exists a constant K depending on m, n, p,q and the 
dimensions of the cone C providing the cone condition for Vt, such that 
for allue W'^'^iQ), 

\\u\\Li{n) < K\\u\\1ym,p(^a)\\''^\\Lv(n) ^ (^-1) 
where 9 = (n/mp) — (n/mq). 

We write K = K{m,n,p, q,Q). We always use Theorem ID. II with 
n = 3, m = 1, and p = p,q = 3p for some p > 3. Hence mp > n, 
p < q < oo, and 6 = 6{p) = 2/p < 1. Moreover, we always use it with 
Q being a ball, whose radius in all cases is bounded from above by 1 
and from below hj R/2 for some R> fixed. 

Let Ko = Ko{p) = K{1, 3,p, 3p, Bi{0)) with Bi{0) C the unit ball 
(which does satisfy the cone condition). Note that then, by scaling, 
fID.ll) implies that for all r < 1 and all xq G M^, 

\H\L^v(Br{^o)) < ^0'""l«Ci.P(B,(xo))ll^llip(V(xo)) ' (^-2) 

with e = 2/p. 

To summarize, we therefore have the following corollary. 

Corollary D.2. Let p > 3 and R G (0, 1]. Then there exists a constant 
K2, depending only on p and R, such that for all r G [R/2,1], xq G M^, 
and allue W^'P{Br{^o)) , 

\\u\\L3p{Br{^o)) < ^2||M||^i,P(i?,(xo))ll^llip(V(xo)) ' (^-3) 

with e = 2/p. 
Here, 

K2 = Mp, R) = {2/Rf'PKo{p) , (D.4) 

where Kq{p) = K{l,3,p,3p, Bi{0)) in Theorem ID. II above. 

Theorem D.3. [H Theorem 4.2] Let Q be a bounded domain in M" 
and let a'^ G C{Tl), b\ce i,j G {!,..., n}, with A,A > such 
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that 

n 

J2 - ^l-^l^ ' xeQ,^eW, (D.5) 

n n 

||a*-'||L'-(Q) + ^ ||&lL°°(n) + ||c||L°°(f7) < A. (D.6) 

i,j=l i=l 

Suppose u G W^^liyt) satisfies 

n n 

Lu=J2- a'^DiDjU + ^ UDiU + cu = f . (D.7) 

i,j=l i=l 

Then for any Q' GG Q, 

\\u\\w^:P(^Q') < C{j\\f\\LP{n) + \\u\\LP{n)} , (D.8) 

where C depends only on n,p,A/X,dist{Q',dQ}, and the modulus of 
continuity of the a*-^ 's. 

We use Theorem ID. 31 in the case where Q' and Q are concentric balls 
(and with n = 3, p = 3p/2, a^^ = 5ij,U = c = 0; hence A = A = 
1). Reading the proof of the theorem above with this case in mind 
(see [4, Lemma 4.1] in particular), one can make the dependence on 
distjfi', dVL} explicit. More precisely, we have the following corollary. 

Corollary D.4. For all p > 1 there exists a constant = K^{p) > 1 
such that for all u E W^'^P/'^iBr+si^o)) (with xq G R^ r, 5 > 

\\u\\w2,3p/2(^Br(xo)) 

< i^3{||A^^||L3p/2(iJ,+,(xo)) +'^''II^IIl3p/2(B,+,(xo))} • (D.9) 

Theorem D.5. [8, Theorem 5, Section 5.6.2 (Morrey's inequality)] 

Let Q be a bounded, open subset in M", n > 2, and suppose dQ is 
. Assume n < p < oo, and u G W^'^{Vt). Then u has a version 
u* G C^'"'{VL), for'-) = 1 — n/p, with the estimate 

ll^*llcO'T(n) - -^4||M||vyi.p(n) • (D.IO) 
The constant depends only on p,n, and Vt. 



l^llc;o.7(n) := sup |m(x)| + sup — - — — — . (D.U) 



= u* 


a.e.. Above 


|m(x) 


-^y)l 


X 


-y|^ 



Of course, sup^-g^ |m(x)| < 



42 A. DALL'ACQUA, S. FOURNAIS, T. 0. S0RENSEN, AND E. STOCKMEYER 



Remark D.6. Note that [8, p. 245] uses a definition of tlie VT^'P-norm 
wliicfi is sliglitly different from ours (see flTK]) ). but whicli is an equiv- 
alent norm by equivalence of norms in finite dimensional vectorspaces. 
Therefore, f ID.lOp holds with our definition of the norm (but the con- 
stant K4 is not the same as the one in ^ Theorem 5, Section 5.6.2]). 
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